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ABSTRACT 

We present the results from Fabry-Perot imaging spectroscopy of the Balmer-dominated supernova 
remnant DEM L 71 (0505—67.9) in the LMC. Spectra extracted from the entire circumference of the 
blast wave reveal the broad and narrow component Ha line emission characteristic of non-radiative 
shocks in partially neutral gas. The new spectra of DEM L 71 include portions of the rim that have 
not been previously observed. We find that the broad component width varies azimuthally along the 
edge of DEM L 71 , ranging from 450±60 kms" 1 along the eastern edge to values as high as 9851^" 
kms -1 along the faint western edge. The latter width is nearly 60% larger than the value determined by 
earlier spectroscopy of the brightest Balmer-dominated filaments. In part of the faint northern rim the 
broad component is not detected, possibly indicating a lower density in these regions and/or a broad 
component width in excess of 1000 kms -1 . Between the limits of zero and full electron-ion temper- 
ature equilibration at the shock front, the allowed range of shock velocities is 430—560 km s" 1 along 
the east rim and 700—1250 kms -1 along other parts of the blast wave. The Ha broad-to-narrow flux 
ratios vary considerably around the remnant, ranging from 0.4 to 0.8. These ratios lie below the val- 
ues predicted by our shock models. We find that narrow component Ha emission from a cosmic ray 
precursor may be the cause of the discrepancy. The least decelerated portions of the blast wave (i.e., 
regions excluding the brightest filaments) are well characterized by Sedov models with a kinetic energy 

^51 = (0.37±0.06) Z?5q 2 , where D50 is the LMC distance in units of 50 kpc. The corresponding age for 
DEM L 71 is (4360±290) I?5o yr. This is the first time that velocity information from the entire blast 
wave has been utilized to study the global kinematics of a non-radiative SNR at a known distance. 

Subject headings: ISM: supernova remnants: individual (0505—67.9, DEM L 71 )-ISM: kinematics and 
dynamics, shock waves 



1. INTRODUCTION 

One distinctive property of non-radiative supernova 
remnants (SNRs) in partially neutral gas is the presence 
of optical line emission which allows us to trace the pre- 
cise positions of their blast waves. The optical spectra 
of these SNRs (termed non-radiative because postshock 
cooling losses are negligible) are dominated by Balmcr line 
emission, produced by collisional excitation when neutral 
hydrogen is overrun by the blast wave (Chevalier & Ray- 
mond 1978, Bychkov & Lebedev 1979). Each emission line 
consists of two components: (1) a narrow velocity com- 
ponent produced when cold, ambient H I overrun by the 
shock is excited by electron and proton collisions, and (2) 
a broad velocity component produced when fast neutrals 
created by postshock charge exchange are collisionally ex- 
cited (Chevalier, Kirshner & Raymond 1980). Since the 
optical emission arises in a very thin (<^ 10 16 cm) ion- 
ization zone, measuring the width of the broad Balmcr 
line yields the proton temperature immediately behind 
the shock. The broad-to-narrow flux ratio, on the other 
hand, is sensitive to both the initial electron-ion tempera- 
ture equilibration at the shock front and the neutral frac- 
tion of the preshock gas (Chevalier, Kirshner & Raymond 
1980, Smith et al. 1991, hereafter SKBW91). Together, 



the broad-to-narrow ratio and broad component width are 
highly useful tools for estimating the shock velocity, degree 
of electron-ion equilibration and even the preshock neutral 
fraction in SNRs (SKBW91, Ghavamian 1999, Ghavamian 
et al. 2001; hereafter GRSH01, Ghavamian et al. 2002; 
hereafter GRHB00). 

One of four known Balmer-dominated SNRs in the 
LMC, DEM L 71 was originally discovered as an unre- 
solved X-ray source by Einstein (Long, Helfand & Gra- 
belsky 1981). The remnant appears as a limb brightened 
shell, tear-drop shaped and approximately 1.4' x 1.2' in 
size. Since its discovery, DEM L 71 has been studied 
in the optical via narrow band imagery (Tuohy et al. 
1982) and low and high resolution longslit spectroscopy 
(SKBW91, Smith et al. 1994, hereafter SRL94). At X- 
ray energies, the remnant has been studied spectroscopi- 
cally with ASCA SIS (Hughes, Hayashi & Koyama 1998; 
hereafter HHK98) and Chandra ACIS-S (Hughes et al. 
2002, hereafter HGRS03; Rakowski et al. 2002, hereafter 
RGH03). In their analysis of the global X-ray spectrum, 
HHK98 assumed a Sedov dynamical evolution and applied 
self-consistent non-equilibrium ionization models to derive 
physical parameters for DEM L 71 and six other remnants 
in the LMC. They found that the Fe abundance in DEM L 
71 was twice as high as that of the other remnants in the 
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sample, leading them to conclude that DEM L 71 is the 
remnant of a Type la explosion. In a more recent Chan- 
dra analysis of spatially resolved ejecta material, HGRS03 
again measured a highly elevated abundance of Fe and 
other Si-group elements, making an even stronger case for 
a Type la origin of DEM L 71 . 

While X-ray analyses have allowed us to build a global 
picture of DEM L 71 , optical spectroscopy has thus far 
been limited to longslit observations passing through the 
bright east/west limbs of the SNR (see Figure 1). In these 
cases the slit was oriented east/west along the diameter 
of the remnant, resulting in detections of broad (^600 
kms -1 FWHM) and narrow Ha emission from the bright- 
est optical filaments and some of the fainter emission be- 
tween them (SKBW91, SRL94). However, both Tuohy 
et al. (1982) and SKBW91 noted that the broad and nar- 
row Ha emission from the bright filaments was accompa- 
nied by moderate [S II], [N II] emission, a sign that the gas 
behind these shocks has begun to cool. In this case, since 
the narrow component Ha emission is produced both by 
collisional excitation at the shock front and recombination 
of the cooling gas, the broad-to-narrow ratios of the bright 
filaments cannot be used to infer the shock parameters. 

Not all of the bright Ha in DEM L 71 is confined to fil- 
aments. Prominent, localized lumps of emission can also 
be seen in the interior. Three particularly bright patches 
around 5"in size are located 10"— 20" inside the west rim. 
(Figure 1). These patches are also visible in narrow band 
[S II] images of DEM L 71 . At least half a dozen other 
lumps can also be seen scattered throughout the interior of 
DEM L 71 . Interestingly, a comparison of our narrowband 
Ha images with Chandra images of DEM L 71 reveals that 
the brightest of the interior patches have no X-ray coun- 
terparts (RGH03). The optical emission from these clouds 
may be produced by clumpy neutral gas which has been 
photoionized by He II A304 emission from the blast wave 
(GRHB00). 

Given the difficulty of separating radiative from non- 
radiative shocks in prior optical spectra of DEM L 71 , 
we undertook a new spectroscopic study to investigate 
the more widely distributed optical emission around this 
remnant. While considerably fainter than the promi- 
nent E/W filaments, the fainter emission extends unbro- 
ken over a larger portion of the circumference and traces 
pure Balmer-dominated shocks, as evidenced by the lack 
of [S II] and [O III] emission from these shocks in the im- 
ages of Tuohy et al. (1982) (and our own unpublished 
imagery). Assuming that variations in Ha surface bright- 
ness reflect density variations in the surrounding ISM, the 
fainter emission should trace portions of the blast wave 
that have suffered the least deceleration. To investigate 
the entire network of Balmer-dominated shocks in DEM 
L 71 , we observed the remnant with the Rutgers Fabry- 
Pcrot imaging spectrometer (RFP), centered on the nar- 
row component Ha line at the velocity of the LMC. Since 
the entire SNR (1.4' x 1.2') lay within the RFP field of 
view, we were able to extract spectra from the entire blast 
wave rim while simultaneously measuring the radius of 
each extraction point from the center of DEM L 71 . This 
was invaluable for correlating the variation of shock speed 
with both radius and position angle along the blast wave 
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Combined with the recent observation of DEM L 71 by 
Chandra (HGRS03, RGH03) the RFP spectra present us 
with an unprecedented opportunity to probe the under- 
lying physics of high Mach number, collisionlcss shocks. 
Thanks to its 0.5" resolution, Chandra has revealed a 
strikingly detailed correlation between the X-ray and op- 
tical morphologies of the blast wave, suggesting a common 
origin for both sources of emission (HGRS03, RGH03). As 
shown by RGH03, the postshock proton temperatures and 
shock velocities measured from Balmer-dominated spec- 
tra can be combined with the postshock electron temper- 
atures of the blast wave measured from Chandra spectra 
to measure the equilibration fraction immediately behind 
the shock ((T e /T p )o) as a function of shock speed. As 
we will show in this paper, the RFP and Chandra data 
probe shock speeds that vary by a factor ~ 3 along the rim 
of DEM L 71 , therefore we can use our multiwavclcngth 
study of this SNR to probe the (poorly known) degree of 
collisionlcss heating over a large sample of shock speeds, 
in a model-independent manner. 

2. SPECTROSCOPIC OBSERVATIONS 

Our observations of DEM L 71 were performed with the 
RFP spectrometer at the 1.5 m telescope of Cerro Tololo 
Inter-American Observatory on 1998 February 18 and 
February 21 (UT). The Tekl024 CCD and 200 mm camera 
lens were used with the f/7.5 secondary focus, giving an 
image scale of 0"65 pixel -1 . The field of view was a cir- 
cle 7.8' in diameter, centered on coordinates 05 h 05 m 42!8, 
-67°52'35"9 (2000). The scans were centered on the LMC 
systemic velocity Vhelio(LMC) = +278 kms -1 and sam- 
pled the Ha line profile at 14 velocity slices A Vi = (—704, 
-617, -480, -347, -210, -155, -128, -77, -18, +87, 
+ 192, +330, +471, +791) kms" 1 . All scans except for 
those at —704 and +791 kms _1 were performed on the 
first night of the observing run. One image was acquired at 
each etalon setting, with an exposure time of 750s at —704 
kms _1 and +791 kms _1 and 500s at all other settings. 
The average seeing during the first night was 2'.'2 FWHM; 
the worst seeing occurred during the scans centered red- 
ward of the Ha narrow component. During the second 
night, the average seeing was 2'.'8 FWHM. In our observa- 
tions we used the CTIO Ha filter 6563/75 (X F whm=75 
A) to isolate emission from a single spectral order. 

We reduced the Fabry-Perot images using IRAF 2 and 
our own custom software. We applied overscan and bias 
subtraction to all images in the standard way. Separate 
flat field images of the telescope dome white spot were 
obtained at each observed wavelength, and were used to 
correct the detector pixel-to-pixel response variations and 
to remove the Ha filter transmission response. The wave- 
length scale was calibrated using a series of comparison 
exposures of Ha and nearby Neon lines. The zero-point 
drift of the wavelength solution was determined from Ha 
calibration exposures interspersed with the object frames 
during the observations. We estimate that the resulting 
wavelength calibration is accurate to better than 0.1 A. 
The RFP instrumental spectral profile measured from the 
calibration exposures is well fit by a Voigt function with 
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a Gaussian width of 44 km s and Lorentzian width of 
52 kms^ 1 . The resulting instrumental FWHM is 117 
kms -1 . We removed cosmic rays from the images by in- 
terpolating over neighboring pixels and used foreground 
stars to register all the images to the one at the central 
velocity, AV = — 18 kms -1 . We scaled the transmission of 
the RFP scans to a common airmass using the CTIO mean 
extinction tables, then applied a Gaussian smoothing ker- 
nel to each image to produce a common effective seeing of 
2'.'8 FWHM. Finally we computed a coordinate solution 
(accurate to within 0'.'2) for each image using USNO— A2.0 
Catalogue stars in the field. 

3. SPECTRAL EXTRACTION AND BACKGROUND 
SUBTRACTION 

The RFP datacube enables us to extract spectra from 
any desired point along the blast wave of DEM L 71 . Us- 
ing snapshots of DEM L 71 at different velocity intervals, 
we can decide on appropriate locations, sizes and orienta- 
tions for spectral extraction apertures. This is particularly 
advantageous in the case of DEM L 71 , where the mix- 
ture of non-radiative and partially radiative shock emis- 
sion around the remnant (Tuohy et al. 1982, SKBW91, 
SRL94) forces us to exercise caution when selecting re- 
gions for spectral analysis. 

In addition to having Ha narrow components contami- 
nated by radiative recombination, most of the bright lumps 
and filaments exhibit [N II] AA 6548, 6583 emission lines 
in their spectra. The broad component widths are large 
enough to overlap with these lines, complicating the task of 
fitting the Ha profiles. Therefore, we focused our spectro- 
scopic analysis on parts of the SNR which have suffered 
less deceleration than those studied earlier by SKBW91 
and SRL94. 

After trying various aperture selection schemes, we de- 
cided that the best option would be to preserve a constant 
number of integrated counts for all of the spectra. This 
results in Ha line profiles with similar signal-to-noise and 
facilitates the inter-comparison of their broad-to-narrow 
ratios. Since the thickness of the blast wave filaments is 
nearly constant throughout the remnant, we chose aper- 
tures with constant widths and varying lengths. While 
selecting spectral extraction apertures for DEM L 71 , we 
were careful to avoid the aforementioned regions showing 
signs of radiative cooling, namely the bright filaments seen 
in the eastern, northern and western portions of the rem- 
nant. We also excluded regions where stars were present. 
The final selection of 14 apertures is marked in Figure 1 
and listed in Table 1. Each aperture contains approxi- 
mately 13,000±500 counts over the entire line profile, be- 
fore sky subtraction. The aperture sizes, locations and po- 
sition angles were chosen to cover as much of the blast wave 
emission as possible while preserving a constant number 
of counts per aperture and avoiding the brighter filaments 
on the east-west sides of DEM L 71 . 

The RFP spectral extraction consisted of summing the 
emission within each defined aperture, for each frame of 
the datacube. Although the wavelength varies with posi- 
tion in each FP image, the variation over each aperture 
is < 0.01 A, negligible compared to the spectral resolu- 
tion. To perform the sky subtraction, we selected a ring 
of pixels in each frame of the datacube, centered on the 
optical axis and passing through the middle of each ex- 



traction aperture. This ensured that the sky fluxes were 
extracted at the same wavelength as the emission within 
each blast wave aperture. The ring included all the sky ly- 
ing off the supernova remnant which did not fall onto stars 
or other detectable discrete sources. We summed the sky 
emission along the ring, rejecting pixels deviating more 
than 4cr from the mean value. We then multiplied the re- 
sulting sky spectrum by N pix {obj) / N pix {sky) to obtain a 
scaled sky spectrum for each of the blast wave apertures. 
Subtracting this spectrum from that of our object yielded 
the final object spectrum for each aperture (Figs. 2a— 2d). 
The sky subtraction removes large scale sources of back- 
ground emission in the RFP data bandpass, namely night 
sky Ha, galactic Ha and [N II] A6583, and the LMC Ha 
and [N II] A6548. 

4. LINE PROFILE FITS 

The intrinsic SNR Balmer line profiles are well described 
by the sum of two Gaussian shapes (Ghavamian 1999) 
with different widths, velocities and strengths. These are 
convolved with the instrumental response function of the 
Fabry-Perot, a Voigt profile. Therefore, the resulting line 
shape is the sum of two Voigt profiles. Since the narrow 
component line is unresolved, it is characterized by two 
parameters: its central wavelength and flux. The broad 
component, on the other hand, is described by three pa- 
rameters: its central wavelength, Gaussian width and flux. 
In fitting the sky-subtracted profiles we leave the five listed 
parameters free, while fixing the baseline level to zero. 

Due to the presence of Lorentzian wings in the extracted 
line profiles, it is not always obvious from the data whether 
the wings of an Ha profile are caused by a second (broad) 
component or whether they are simply the Lorentzian tail 
of a single, narrow component. Since the broad-to-narrow 
flux ratio can vary significantly from point to point along 
the rim of DEM L 71 , we required a method for quan- 
titatively measuring the likelihood of a broad component 
detection in each extracted spectrum. We used the fol- 
lowing approach. First, we fit each Ha profile with a sin- 
gle, narrow Voigt function. Using parameters output from 
a ^-minimized fit, we then added a second Voigt func- 
tion and re-computed the best fit. The initial guesses for 
the broad component parameters were to take Ib/In = 
0.5 and to assume an unresolved broad component profile, 
a G (b) = <J G (n). 

To statistically evaluate the influence of a second Voigt 
component on the profile fits, we computed an F-statistic 
for each profile using the definition of F and the goodness- 
of-fit parameter P(> F) from Martin (1971) and Band 
et al. (1996). In this procedure a large value of F and small 
value of P(> F) would indicate that it is unlikely the supe- 
rior two-component fit resulted from a random statistical 
fluctuation. As Table 2 shows, P(> F) is on the order of 
5% for all blast wave apertures other than apertures 1 and 
2, indicating that a second, broad component is detected 
in each Ha profile. 

5. BLAST WAVE SPECTRA 

Most of the blast wave around DEM L 71 follows a 
'faceted' morphology, where relatively straight sections of 
the shock front appear joined together at abrupt breaks 
in orientation angle. This is most strikingly seen in the 
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northern, southwestern and southeastern edges of the rem- 
nant. The size of each section reflects the length scale for 
variations in ambient density around DEM L 71 . 

The Ha profile fits in Figures 2a— 2d show that the 
Ha broad component is detected in 14 of the 16 spectra. 
The fitting results are summarized in Table 2. The broad 
component widths and broad-to-narrow ratios vary con- 
siderably with position along the rim. The faintest blast 
wave filaments tend to exhibit the largest broad compo- 
nent widths (with the notable exception of aperture 10), 
while sections of the blast wave lying adjacent to filaments 
of intermediate brightness (such as apertures 4, 13 and 14) 
tend to exhibit smaller broad component widths. 

Spectra from apertures 1 and 2 at the northern rim are 
well fit by a single narrow Ha emission line, indicating 
that the broad component emission either lies beneath 
the detection threshold or is entirely absent at these lo- 
cations. Interestingly a broad component is detected in 
the aperture 3 spectrum, even though it is part of the 
same straight, faint filament as apertures 1 and 2 (Fig- 
ure 1). This is likely due to a gradient in preshock neu- 
tral density which causes the surface brightness to increase 
by nearly 70% from aperture 1 to aperture 3. A density 
gradient is also apparent in the X-ray surface brightness 
of the Chandra data (RGH03), indicating that the total 
preshock density also rises between apertures 1 and 3. The 
broad component widths of apertures 3 and 4 and their 
corresponding errorbars are consistent with a single broad 
component width of 840^ioo kms -1 , nearly 40% larger 
than the widths measured by SKBW91 from the brighter 
east/west filaments. 

Progressing along the western side of DEM L 71 , we find 
that the spectra of apertures 5 and 6 exhibit the largest 
broad component widths (~ 1000 kms -1 ) seen to date in 
this SNR. The faint blast wave bulges out beyond the 
brightest western filament (Figure 1). Since the Balmer 
emission is produced very close to the blast wave, the ra- 
dial separation between two filaments is entirely due to 
geometric projection. This suggests that the brightest 
western emission marks the location where the forward 
shock has slowed due to an encounter with denser mate- 
rial. This is confirmed by both longslit spectra (SKBW91) 
and our RFP profiles, which indicate that the broad com- 
ponent width of the bright western filament is only half 
that of apertures 5 and 6. The broad component profiles 
of apertures 5 and 6 are consistent with a single FWHM 
of 985^0 kms- 1 and I B /I N = 0.54±0.09. These are the 
largest broad component widths reported so far in DEM 
L 71. 

The southern edge of DEM L 71 reveals some of the 
best examples of the 'faceted' rim morphology mentioned 
above. A noticeable brightening occurs in the filament 
covered by apertures 7—9. The Chandra image of DEM 
L 71 shows that the enhanced Ha is accompanied by a 
sharp rise in X-ray emission (RGH03, HGRS03). This 
indicates a higher total preshock density at the regions 
covered by apertures 7—9. All three apertures show sim- 
ilar broad component widths and broad-to-narrow ratios 
(Table 2). Combining the spectra from the three aper- 
tures, we obtain a broad FWHM of 805^° kms" 1 and 
Ib/In = 0-49^Q;Qg, similar to values seen in apertures 3 
and 4. The broad-to-narrow ratio and broad component 



width vary considerably from aperture 9 to 12, suggesting 
substantial variations in shock velocity, preshock density 
and possibly preshock neutral fraction. The component 
widths of the aperture 13 and 14 spectra are nearly equal 
to the values reported by SKBW91 for the bright, partially 
radiative eastern filament. 

Unlike the other shocks in our study, the aperture 15 
and 16 filaments are seen well inside the blast wave rim 
(Figure 1). The filaments sampled by these apertures 
are of similar surface brightness, yet the broad compo- 
nent width and broad-to-narrow ratio of the aperture 15 
spectrum are more than twice as large as the aperture 16 
values. The broad component shift in aperture 16 is 160 
km s _1 relative to the narrow component, suggesting that 
the Ha emission in this aperture is dominated by partially 
face-on shocks. The Chandra image of DEM L 71 shows 
that there are in fact three nested blast wave shocks in 
the southwest region (RGH03). By comparing the optical 
images with the Chandra image we have determined that 
two of these nested shocks are sampled by aperture 15. In 
that case, variations in the line-of-sight bulk velocity from 
inter-filament emission could contribute substantially to 
the broad component width of the aperture 15 spectrum. 

6. SHOCK MODELS 

To model the Balmer line profiles, we used the 
one-dimensional, plane-parallel shock code described in 
Ghavamian (1999) and GRSH01 to compute a grid of 
numerical models. The code calculates the density and 
temperature of electrons, protons and hot neutrals behind 
a Balmer-dominated shock. Broad-to-narrow flux ratios 
are computed by a Monte Carlo simulation which follows 
the excitation of Ha photons in the broad and narrow 
components. The radiative conversion of Ly (3 photons 
into Ha (Ly /3 trapping) is also followed in the simula- 
tion. At the shock speeds relevant to DEM L 71 (<; 1500 
kms -1 ) the proton— H ionization and excitation rates are 
< 1% of the electron rates. Therefore, all direct collisional 
ionization/excitation processes are dominated by the elec- 
trons. The shock velocity, fractional electron-ion equili- 
bration (T e /T p )o at the shock front by plasma turbulence, 
and preshock neutral fraction of H are taken as free pa- 
rameters. The quantity (T e /T p ) varies from m e /m p (~ 0) 
for zero equilibration at the shock front to (T e /T p )o = 1 
for full, prompt equilibration at the shock front. (In this 
work we have chosen to describe our results in terms of 
(T e /T p ) rather than the related quantity f eq used in ear- 
lier works (GRSH01, Ghavamian et al. 2002), to main- 
tain a consistent notation between this paper and the com- 
panion work by RGH03). 

In addition to the shock code improvements listed in 
Ghavamian ct al. (2002), we have included the energy 
lost by electrons to collisional excitation and ionization. 
Loss terms are included for collisional ionization of H and 
He, and for Lya and two-photon continuum excitation of 
H and He. Clearly, ionization losses are most important at 
the lowest shock speeds, lowest equilibrations and highest 
preshock neutral fractions, where the electron temperature 
close to the shock front is lowest. In these cases T e can 
initially decline with position behind the shock as excita- 
tion and ionization losses exceed energy gain by Coulomb 
collisions. Subsequently as the ionization fraction rises, T e 
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also rises over a length scale comparable to the thickness of 
the H ionization layer as Coulomb collisions come to dom- 
inate the energy balance. In models with shock speeds 
~500 kms -1 , (T e /T p )o <^0.4 and preshock neutral frac- 
tions J> 0.9, we estimate that ionization losses produce a 
5%— 7% increase in the broad-to-narrow ratios relative to 
models that neglect these losses. 

Using the method outlined in GRSH01, we computed a 
grid of shock models for spectra extracted from the blast 
wave rim. We matched each grid point in (T e /T p )o with 
the shock velocity Vs required to match the observed 
FWHM of the broad component line. Using this set of 
parameters (Vs , (T e /T p )o), we calculated Ib/In^ot a 
range of assumed preshock neutral fractions /#o , then 
compared the results with the observed I B /I N . How- 
ever, before producing the grid we decided to combine 
spectra of regions exhibiting consistent broad-to-narrow 
ratios and broad component widths (to within the errors). 
The main purpose of combining the spectra was to de- 
fine large enough apertures along the blast wave to ob- 
tain adequate statistics in our parallel analysis of Chandra 
data (RGH03). The five combined regions are labeled 3+4, 
5+6, 7+8+9, 11+12 and 13+14 in reference to the origi- 
nal 14 apertures (c.f. Table 2). The combined apertures 
correspond to X-ray spectral extraction regions XI— X5 
from RGH03. After combining the spectra, we re-fit the 
line profiles in the manner described earlier to estimate 
the broad component FWHM and broad-to-narrow ratios 
(Tabic 3). The limits on Vs implied by the combined aper- 
ture broad component widths are listed in Table 4 for the 
cases (T e /T p ) = m e /m p and (T e /T p ) =1. 

7. COMPARISON BETWEEN OBSERVATIONS AND SHOCK 
MODEL PREDICTIONS 

From Figure 3 it is obvious that there is a major dis- 
crepancy between the observed and predicted broad-to- 
narrow ratios: the observed ratios lie as much as 50% be- 
low the smallest modeled values. Clearly, some aspect of 
the shock physics is not included in the models. It seems 
likely that our models overpredict the broad-to-narrow 
ratios because they undcrprcdict the narrow component 
Ha emission. Narrow component Ha emission can be en- 
hanced by processes not included in our Balmcr-dominated 
shock models. Although our calculations include the en- 
hancement of narrow component Ha from Ly/3 trapping 
within the shock, other possible sources of narrow Ly/3, 
such as other shocks in the SNR, are ignored. In addi- 
tion, if the preshock gas is heated in a precursor it may 
produce enough collisionally excited Ha to lower the ob- 
served broad-to-narrow ratios. 

One way of bringing the model broad-to-narrow ratios 
into agreement with the observations would be to forgo 
the Monte Carlo simulation altogether and simply com- 
pute the integrated Ha flux assuming a negligible Ly/3 op- 
tical depth in the broad component (case A conditions) 
and a high Ly/3 optical depth in the narrow component 
(case B conditions) (Chevalier, Kirshner & Raymond 1980, 
SKBW91, Ghavamian 1999). However, we find that al- 
though the Ly/3 optical depth can be large enough at line 
center to make case B a valid approximation for the narrow 
component, it is generally not small enough at broad com- 
ponent line center to allow a case A approximation for the 
broad component. For example, tl v 8 (narrow) <~ 7 and 



T Lvf3 (broad) w 0.2 for V s ~ 1000 kms _1 and f H o= 0.9. 
Smaller neutral fractions can lower TL y [j(broad), but not 
enough to ignore Ly/3 trapping in the broad component. 
Therefore, usage of the case A/case B approximations in 
the computation of broad-to-narrow ratios requires physi- 
cal conditions not met in practice. 

The limiting assumptions from above can be made more 
plausible if we include absorption of narrow Ly/3 photons 
produced in other parts of the remnant. For a given shock, 
this would effectively increase TL y /3(narrow). Absorption 
of nonlocally produced Ly/3 is more important for a rem- 
nant like DEM L 71 which is completely surrounded by 
Balmer filaments than remnants with only partial shells 
of Balmer emission. A Ly/3 photon may propagate all the 
way to the far side of the remnant without being absorbed. 
Once the photon enters the shock structure on the far side 
of the remnant, it will most likely be converted into nar- 
row Ha as it propagates upstream. From the Monte Carlo 
shock models we estimate that for shock speeds between 
600 and 1000 kms -1 , preshock neutral fractions greater 
than 0.1 and all values of (T e /T p ) , less than 8% of all 
narrow component Ly/3 photons excited behind the DEM 
L 71 blast wave escape downstream. If all of these pho- 
tons are absorbed in other parts of the blast wave and are 
converted into narrow Ha, the broad-to-narrow ratios be 
lowered by at most 20% from their intrinsic values. This 
is less than half the amount required to match most of the 
observed broad-to-narrow ratios with the smallest model 
values. Although the Ly/3 luminosity of the partially ra- 
diative shocks in DEM L 71 is nearly twice as large as that 
of the fainter, pure non-radiative shocks, their narrow Ly/3 
flux is offset by the significantly smaller surface area of the 
remnant covered by the partially radiative shocks. There- 
fore, we now consider the possibility that most of the ex- 
cess narrow component Ha is produced ahead of the blast 
wave, in a precursor. 

8. EXPLAINING THE BROAD-TO-NARROW RATIOS 

When ionizing photons or energetic particles produced 
behind a non-radiative shock (such as cosmic rays or fast 
neutrals) cross upstream, they generate a layer of heated 
preshock gas known as a precursor. The characteristic 
size d p of the precursor depends on the mean free path 
traveled by the photon or particle before depositing its 
energy into the upstream gas. This energy is shared be- 
tween electrons and ions via Coulomb collisions, which for 
the preshock temperatures inferred from spectroscopic ob- 
servations (12,000 K-40,000 K, Hester, Raymond & Blair 
1994, SRL94, GRHB00) leads to T e = T t ahead of the 
shock. Collisional excitation by electrons within the pre- 
cursor can produce observable line emission (Fesen & Itoh 
1985, GRHB00). The detectability of such emission de- 
pends on the thickness, temperature and ionization struc- 
ture of the precursor. Therefore, we must constrain these 
parameters in order to investigate the influence of precur- 
sor emission on the observed broad-to- narrow ratios. 

A constraint on the precursor thickness can be obtained 
from the optical images of DEM L 71 . If the precursor 
is spatially resolved, it should appear as a layer of diffuse 
Ha emission extending ahead of the Balmer-dominated fil- 
aments (GRHB00). Since the shock front continually over- 
runs the heated ambient gas, the diffuse precursor emission 
should drop sharply behind the Balmer filaments. A nar- 
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row band image of DEM L 71 acquired before the RFP 
observations does show marginal evidence for diffuse emis- 
sion above the local background, extending a few arcsec- 
onds ahead of the northern limb. However, elsewhere in 
the remnant the surrounding sky does not show any ob- 
vious diffuse emission above the background LMC level. 
There is also no clear evidence of diffuse emission ahead 
of most of the filaments in the RFP scans, aside from the 
emission mentioned earlier ahead of the northern rim. 

The remaining possibility is that the precursor emission 
is produced on small scales, comparable to or less than 
the Balmer filament thickness. One method of constrain- 
ing the thickness of such a layer is to examine the RFP 
images for radial stratification between shocks observed 
at pure broad component velocities and shocks observed 
at narrow component velocities (e.g., less than the instru- 
mental width of 117 kms -1 ). If pure broad component 
filaments appear inside of pure narrow component fila- 
ments, this would be nominal evidence for emission ahead 
of the blast wave on scales slightly larger than the Balmer 
filament thickness. Performing this test with the RFP im- 
ages, we find no obvious spatial shift between the narrow 
and broad emission. Therefore, the layer of excess narrow 

component emission is restricted to a thickness <; 3 (the 
seeing FWHM). This restriction is important because it 
allows us to constrain the type of precursor producing the 
excess Ha emission. 

Before examining the various precursor mechanisms in 
detail, we must consider the relationship between the spa- 
tial resolution of the detector and the relative contribu- 
tions of the shock and precursor emission from the Balmer- 
dominatcd filaments. At a given pixel scale, the spatial res- 
olution of a detector imaging an LMC remnant is over an 
order of magnitude worse than that of a detector imaging 
a Galactic remnant. This can make the spatial separation 
of the precursor and shock emission much more difficult 
for an LMC remnant like DEM L 71 . If the thickness of 
the precursor d p and the thickness of the postshock H ion- 
ization zone d s are both smaller than the projected size of 
the seeing FWHM, it will be nearly impossible to subtract 
the precursor emission from that of the shock in a Balmer- 
dominatcd filament. Furthermore, the ratio of precursor 
to shock surface brightness I p /I s depends on the viewing 
angle to the shock. If d p > d s , but both are smaller than 
the projected size of the seeing FWHM, I p /I s can become 
significant for edge-on filaments, even if the ratio is small 
at face-on positions. 

Assuming a spherical geometry and that the radius R of 
the remnant extends from the remnant center to the blast 
wave shock front, the edge-on emission is related to the 
face-on emission through a factor sec(9) = L/d, where 
L is the depth of emission viewed edge-on and d is the 
larger of the layer thickness and the projected size of the 
atmospheric seeing. If d s and d p are both smaller than the 
projected seeing, the ratio is 



(i) 

where we have used R » d s , d p to simplify the expression. 
The ratio of precursor to shock emission for the face-on 



position is 

/Ip\ = n e g Ha (T e )d p 

\I.Jf ~ 0.2V S " 1 ] 

where is the Ha collisional excitation rate in the pre- 
cursor, T e is the electron temperature in the precursor, 0.2 
is the number of Ha photons excited per H atom, Vg is 
the shock speed and n e = fn+no is the preshock electron 
density. Using Equations 1 and 2 and our limits on the 
shock parameters, we can first calculate (I p /I s )f, then use 
Equation 1 to find (I p /I s ) e , the ratio expected for the fil- 
aments. Although some spatial variation of T e and f H + 
is expected within the precursor (Boularcs & Cox 1988, 
GRHB00), we take these quantities to be constant for the 
crude estimates in this section. 

The amount of precursor emission required to match the 
observed broad-to-narrow ratios depends on the assumed 
shock parameters. For an observed broad-to-narrow ra- 
tio (Ib/In)o and specified model broad-to- narrow ratio 
(Ib/In )™, the implied ratio of precursor to shock emis- 
sion is 

Ip = ({lB/lN) m /{I B /lN)o) - 1 , , 

I, " (Ib/InU + 1 [ ' 

where we have assumed that the observed narrow com- 
ponent flux is given by the sum of fluxes from the shock 
and the precursor. From Figure 3, the smallest broad- 
to-narrow ratios predicted by our models lie in the range 
0.9—1.1. The observed broad-to-narrow ratios lie in the 
range 0.4—0.75. Therefore, using Equation 3, the edge-on 
ratio of precursor to shock emission must lie in the range 
0.1 < (I p /I s )e < 0.9 to resolve the broad-to-narrow ratio 
discrepancy. Therefore, matching the smallest (Ib/In)o 
with (Ib/In )m can require nearly equal contributions to 
the Ha flux from the precursor and shock. 

For the purpose of obtaining order-of-magnitude esti- 
mates of Ip/I s , the main distinction between the precur- 
sor scenarios lies in the estimated values of d p and T e . We 
consider three scenarios below: 

(1) Photoionization Precursor— The gas behind all non- 
radiative shocks produces ionizing radiation which propa- 
gates far upstream and heats the preshock gas (Draine & 
McKee 1993, SRL94, GRHB00). The dominant source of 
ionizing photons is He II A304 (Ly a) line emission, which 
can heat the preshock gas enough to produce observable 
optical emission (T e - 12,000 K-20,000 K, GRHB00) but 
not enough to establish ionization equilibrium. The char- 
acteristic size of the photoionization precursor is one mean 
free path of a He II A304 photon, d 3n4 = (n H o o^) -1 
(SRL94, GRHB00). Assuming a total preshock density 
no <^ 0.5 cm~ 3 (as constrained by Chandra spectral fits of 
RGH03) and neutral fractions lower than 0.9, we obtain 
an angular size of > 7" for the precursor, sufficiently large 
to be resolved in the RFP images. However, as mentioned 
earlier, there is no clear evidence of such emission ahead 
of the Balmer filaments and above the LMC background. 
Therefore, we conclude that although a He II A304 is un- 
doubtedly present, it does not contribute enough to the 
Ha filament emission to account for the anomalously low 
broad-to-narrow ratios. 

(2) Fast Neutral Precursor— Fast neutrals on the tail 
of the distribution can cross upstream and deposit en- 
ergy into the preshock medium via clastic collisions and 
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charge exchange with the preshock ions (Hester, Raymond 
& Blair 1994, SRL94, Lim & Raga 1996). Although there 
is considerable uncertainty regarding the energy transfer, 
existing studies suggest that heating in the precursor can 
be substantial (SRL94, Lim & Raga 1996), with tempera- 
tures in excess of 20,000 K easily attained. However, the 
thickness of the fast neutral precursor is on the order of the 
charge exchange mean free path, d p = d cx ~ 10 15 /n « 
2xl0 15 cm for the deduced preshock parameters. We 
estimate the electron temperature in the precursor from 
measured narrow component widths of SRL94. Assuming 
purely thermal line broadening and (T e /T p )o = 1 in the 
preshock gas, the 30 kms -1 — 45 kms -1 narrow compo- 
nent widths quoted by SRL94 imply an electron temper- 
ature of 20,000 K-40,000 K for the precursor. The colli- 
sional excitation rate coefficient qn a rises steeply with T e 
over this temperature range. From the Case B calculation 
of Aggarwal (1983), we estimate q Ha w (0.1-3.5) x 10" 10 
cm~ 3 s^ 1 . We adopt no ~0.5 cm -3 and take the preshock 
ionization fraction f H + to lie between 0.2 and 0.9. Pho- 
toionization by He II A304 and X-ray photons from the re- 
verse shock would make ionization fractions smaller than 
0.2 unlikely, while fractions larger than 0.9 would make 
the Balmer filaments too faint to observe. We assume 
Vs = 900 kms -1 using the average speed listed in Table 4 
for the least decelerated portions of the blast wave. 

Inserting the estimates from above into Equation 2, 
and assuming that d p w d cx , we obtain (I p /I s )f <^0.1. 
Therefore, explaining the full range of broad-to-narrow 
ratios measured in the filaments remnant would require 
{dp/ds) 1 / 2 ;> 10. However, for the shock speeds found in 
DEM L 71 , d s <~ d cx ~ d p , so there is little enhancement 
in the ratio of precursor to shock emission at edge-on view- 
ing angles, (I p /I s ) e « (I p /I s )f. Therefore, narrow Ha 
emission from a fast neutral precursor might contribute to 
the blast wave spectra of DEM L 71 , but not enough to 
resolve the discrepancy in all the broad-to-narrow ratios. 

(3) Cosmic Ray Precursor— Cosmic rays are accelerated 
by the shock when they scatter from self-generated Alfvcn 
turbulence upstream. Since the size of the cosmic ray pre- 
cursor depends on the poorly known diffusion coefficient 
of cosmic ray ions, we can only place an upper limit on its 
size by requiring that the precursor layer be thin enough 
to avoid ionization equilibrium (otherwise no neutral H 
would remain to produce the Balmer-dominated emission). 
This requires d p < vs/n e qi 1 where qi is the ionization 
rate coefficient in the precursor (Hester, Raymond & Blair 
1994, SRL94). Using & (40, 000 K) «3xl0- 10 cm 3 s" 1 
(Janev et al. 1987), we obtain d p <6xl0 17 // ff +. For 
f H + ;>0.2, d p is smaller than the projected size of the 
seeing FWHM. Inserting this upper limit d p into Equa- 
tion 2 and using d s « 10 15 /n H + (the thickness of the H 
ionization zone behind the shock), we find that (I p /I s )f 
can reach values as large as 5.8. According to Equation 1, 
this ratio is enhanced by another factor (dp/dg) 1 / 2 ^25 
at edge-on viewing angles. Therefore, even if d p were sev- 
eral orders of magnitude smaller than the value set by the 
ionization length argument, the narrow Ha emission from 
a cosmic ray precursor would still be able to lower the in- 
trinsic broad-to-narrow ratios down to the observed levels. 

If a cosmic ray precursor exists in DEM L 71 , it should 
generate optical forbidden line emission. Using the pre- 
cursor density, temperature and size from above, we can 



place an upper limit on the [S II] and [N II] surface bright- 
ness. If the preshock ionization fraction is 0.2 <J f H + <J 
0.9, and f s + =1 and log(S/H) + 12 = 6.87 (Russell & 
Dopita 1990), then using the collision strengths of Cai & 
Pradhan (1993) and assuming the low density limit, we 
find that the ratio of precursor [S II] (6716 + 6731) to 
precursor Ha is <; 5xl0~ 3 , too faint for detection. Simi- 
larly, if we take fg+ « f H + and use log(N/H) + 12 = 6.55 
(Russell & Dopita 1990), the ratio of [N II] (A6548+6583) 
to Ha is < 3xl0~ 4 in the precursor. At a temperature of 
40,000 K, the H collisional excitation rate is more than an 
order of magnitude larger than the rate at lower tempera- 
tures (12,000 K— 15,000 K, as found in a photoionization 
precursor, GRHB00). On the other hand, the forbidden 
line excitation rates do not increase as rapidly between 
12,000 K and 40,000 K. Therefore, the Balmer' line emis- 
sion dominates over forbidden line emission in the precur- 
sor. This is consistent with the lack of detection of [S II] 
emission from the pure Balmer-dominated filaments in our 
narrow band images of DEM L 71 . 

Given the uncertain contribution of precursor emission 
to the narrow component in DEM L 71 , we refrain from 
further attempts to use the observed broad-to-narrow ra- 
tios to estimate (T e /T p ) and Vs . Instead, we have com- 
bined the proton temperatures and range of shock speeds 
implied by broad component widths around the rim with 
electron temperature measurements from Chandra data to 
estimate (T e /T p ) and V s (RGH03). In the rest of this pa- 
per, we use these estimates to draw conclusions regarding 
the age and explosion energy of DEM L 71 . 

9. EVOLUTIONARY MODELS OF DEM L 71 

The large radius (~ 10 pc) and moderate blast wave 
speed (~400-1000 kms" 1 ) of DEM L 71 suggest that 
this remnant is well into the Sedov stage of evolution. One 
way to test this conclusion is to compare scale quantities 
predicted by Sedov models, such as the characteristic ra- 
dius, with the observed radius of the supernova remnant. 
In the case of a SNR in a uniform density medium and an 
exponential ejecta profile (as expected for Type la explo- 
sions; see Dwarkadas & Chevalier 1998, Wang & Chevalier 
2001), the blast wave initially evolves through an ejecta- 
dominated phase. It then enters the Sedov stage when its 
radius has grown to a characteristic size Rst, given by 

(Wang & Chevalier 2001), where M e j is the ejected mass, 
M c h is the Chandrasekhar mass (1.4 M@) and no is the 
density of the ambient medium. Taking a range 0.3 <J 
n <;l-4 from fits to the combined aperture Chandra blast 
wave spectra (RGH03; also see Table 4) and assuming 
M e j = M c h, we find Rst ^3.3 pc. This is less than 
half the average blast wave radius of DEM L 71 ; therefore 
the remnant should be well into the Sedov- Taylor phase of 
evolution. 

Two fundamental parameters of interest in DEM L 71 
are the age r and explosion kinetic energy E$i. We as- 
sume that azimuthal variation in Vs and R result from 
azimuthal variations in ambient density. A reasonable first 
approximation for modeling this dependence is to describe 



8 



the blast wave evolution in each spectral extraction aper- 
ture in terms of its own Sedov solution. In that case, E 5 i 
and r are identical everywhere along the blast wave rim, 
with n varying from section to section. In that case, the 
equations governing the evolution of the remnant predict 
Vs = aR, where a = g^. Furthermore, inverting the 
Sedov- Taylor relation (Sedov 1959) for R(t) yields 

n (R) = 3.1 x 1CT 3 t 2 £ 5 i iT 5 cm" 3 , (5) 

where no is the preshock H density and r is measured in 
years. Here, the total preshock density in the Sedov equa- 
tion is related to the preshock H density (measured from 
the X-ray spectra) by a factor of 1.4. Fitting Vs(R) and 
no(R) to the data yields r and E$\. To perform these fits, 
we utilized the values of (T e /T p ) determined by RGH03 to 
calculate the appropriate shock speeds consistent with the 
broad component widths found in the combined aperture 
regions (see Table 4). We also used n (R) determined 
from the X-ray models of RGH03. Since GRSH01 and 
RGH03 found evidence that (T e /T p ) and Vs are inversely 
related, the faster shocks in DEM L 71 will have shock 
speeds closer to the (T e /T p )o = predictions, while the 
slower shocks will have speeds closer to the (T e /T p )o = 1 
predictions. The assumption that the combined aperture 
shocks are either all unequilibrated or fully equilibrated is 
not strictly valid; however, the results of our fits to Vs(R) 
and no(R) do not change significantly if (T e /T p )ois set to 
cither extreme for all the shocks. 

The analysis outlined above requires an estimate of R 
for each aperture, which in turn requires an estimate of 
the expansion center. Therefore, the following procedure 
was implemented: first we used a range of expansion cen- 
ter coordinates to calculate a grid of fits for Vs(R). After 
choosing the expansion center coordinates which yielded 
the smallest % 2 in the grid, we then computed the lcr er- 
ror box about that best fit center (see Figure 1). Next, we 
independently performed the same procedure for no(R), 
obtaining an independent estimate of the expansion center 
and a la error box on that fit (also marked on Figure 1). 
We chose the central pixel from the overlap region between 
the Vs(R) and no(R) error boxes and designated the co- 
ordinates of that pixel (x cen , y ce n) to be the expansion 
center of the remnant. Finally, returning to the grids in 
Vs(R) and no(R), we use the fits at the expansion center 
{xcen, Veen) to estimate r and E 51 . 

In Figure 4 we show the fitted blast wave speeds vs. 
radii, derived using our best estimate for the expan- 
sion center (located at a(2000) = 05 h 05 m 42?3, <5(2000) = 
— 67°52'41'/3). Here we have ignored the uncertainty in 
blast wave radius in our images. Since the distance to 
the LMC is uncertain by as much as 13% (Gibson 2000), 
we have included a factor D 50 (the distance in units of 
50 kpc) when quoting the fitted parameters. The slope 
a of the fit is (89.6±5.9) D^ kms _1 pc _1 , which yields 
r = (4360±290) D 50 years for DEM L 71 . This age is less 
than half the value quoted by SKBW91, whose analysis 
was based on broad component width measurements of the 
slower, brighter Balmcr-dominated shocks. We also note 
that an age of (4360+290) D 50 years is in good agreement 
with the age of 4700 years calculated from fits to ASCA 
data by HHK98. 

Using an expansion center fixed at (x cen , y cen ) and 



our estimated remnant age, the fit to no(R) (shown 
in Figure 5) yields an explosion kinetic energy of 

5 /2 

E51 = (0.37±0.06) D 5 q . The two main sources of uncer- 
tainty lie in the distance to the LMC and the FWHM mea- 
surement of the broad component width (which affects r 
through Vg ) . Allowing for the distance uncertainty to the 
LMC and the statistical error, E 5 i lies between 0.2 and 
0.6, less than half the value computed by HHK98. The 
main reasons for the difference arc that HHK98 approxi- 
mated DEM L 71 as a sphere 10.4 pc in radius (despite the 
sizable variation of R with azimuth) and fit the global X- 
ray spectrum of both the blast wave and metal rich ejecta 
to generate Sedov models for the remnant. Our optical 
study utilizes velocities, densities and radii from multiple 
sections of the blast wave; therefore, it provides a much 
more direct and accurate estimate of the explosion en- 
ergy than the models of HHK98. However, our derived 
explosion energy for DEM L 71 is low compared to val- 
ues predicted by 1-dimensional (Hoflich & Khokhlov 1996) 
and 3-dimensional (Reinecke et al. 2002) hydrodynami- 
cal models of Type la SNe. The reason for the difference 
is not clear; one possible explanation is that some of the 
kinetic energy imparted to the ejecta in DEM L 71 was 
lost during an earlier epoch when, unlike now, cosmic ray 
acceleration was efficient. 

The expansion center overlap region in Figure 1 is 
shifted slightly to the east of the apparent center. This 
is consistent with a supernova explosion in a medium with 
an azimuthal density gradient. A higher than average den- 
sity has decelerated the blast wave on the entire east side 
of the remnant, while recent encounters with H I clouds 
on the east and west sides have caused a strong bright- 
ening in the optical emission. The presence of forbidden 
line emission in these filaments indicates that these shocks 
are becoming radiative. In contrast, the blast wave is ex- 
panding most rapidly along the north and southwest rims 
of the remnant, where the density is lower than average. 
The fact that the preshock densities follow an R -5 rela- 
tion so well implies that despite the large differences in 
shock speed between the east rim of DEM L 71 and the 
rest of the remnant, the kinematics of the blast wave are 
consistent with Sedov evolution. 

One interesting property of DEM L 71 is that although 
the partially radiative filaments on the east rim are nearly 
twice as bright as the filaments covered by apertures 13 
and 14, they all lie at the nearly the same radius and 
have similar broad component widths (as verified by our 
own data and that of SKBW91). The preshock gas in 
these partially radiative, [N II]- and [S II]-cmitting clouds 
is nearly twice as dense as that of the shocks in apertures 
11+12 and 13+14. If the densities of these clouds were 
plotted on Figure 5, they would lie well above the fitted 
curve for no(R). This would reflect a significant departure 
from the Sedov description in the brightest filaments, and 
would suggest that while the shocks in apertures 11+12 
and 13+14 are moving through a moderate, large scale 
density gradient, the brightest filaments are shocks driven 
recently by the blast wave into stronger, more discontinu- 
ous jumps in ambient density. 

In the above analysis we have assumed that the blast 
wave contained within each combined aperture has propa- 
gated through a constant density ISM during the course of 
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the remnant's evolution. This is a reasonable assumption 
since the accumulated weight of evidence suggests that 
DEM L 71 was produced by an exploding white dwarf 
(Type la SN). However, we can also determine the age 
and explosion energy of DEM L 71 under the assumption 
that the remnant was produced by a Type lb or Type II 
SN inside of a low density cavity. In that case, the bubble 
would have been evacuated by a strong wind from a mas- 
sive (;>8 M Q ) prognitor star. The density profile of the 
cavity would follow a power law p(R) oc , where u) < 
(Cox & Franco 1981, Franco et al. 1991). Taking u = -2 
(Franco ct al. 1991), the relationship between Vs and R 
becomes Vs = % — ■ Using Equation (14) of Cox & Franco 
(1981), the coefficient of R~ 5 in Equation 5 is multiplied 
by (1 - 0.239w)/(3 - w). The R~ 5 proportionality in 
Equation 5 remains unchanged because the R 2 dependence 
of Uq{R) is already contained in our measurement of the 
current density from the X-ray data. 

Taking the fitted slopes for Vs(R) and no(R) quoted 
in the previous subsection, we estimate that for the case 
lu = -2 the age of DEM L 71 is r = (3100±200)D 50 yr 

and £ 5 i = (0.33±0.05)D 5 5 /2 . Thus both the implied ex- 
plosion energy and the estimated age of the remnant are 
slightly lower when the blast wave propagates through a 
cavity. This result is understandable because less energy 
and less time is required to expand the SNR to its cur- 
rent radius if the blast wave has only recently encoun- 
tered the walls of the cavity. Although steepening the 
density profile to u) < — 2 would lower the blast wave en- 
ergy still further, an asymptotic limit is quickly reached 

5 /2 

where E 51 ->0.27D 50 ' as lo — > — 00. 

The assumption that each section of the DEM L 71 blast 
wave follows its own Sedov evolution may be questioned on 
the grounds that pressure gradients transverse to the blast 
wave could cause one section to affect the evolution of an- 
other section. In the general case, the blast wave should 
be modeled using a multidimensional hydrodynamic de- 
scription. However, as indicated in Table 4, only the 
shock parameters and radii from region 13+14 differ signif- 
icantly from those of the remaining regions. Therefore, the 
strongest pressure differences should arise between regions 
on opposite sides of DEM L 71 . The relevant question, 
then, is to determine how the time T cross required for a 
sound wave to travel across the remnant interior compares 
to the age r of the remnant. Using the self-similar rela- 
tions of Cox & Franco (1981) to compute the sound travel 
time as a function of average interior temperature, we then 
integrate the relation over radius to obtain T cross . Taking 
the current blast wave speed to be 1000 kms" 1 and the 
diameter of the remnant to be 20 pc, we find that for a 
remnant age ~ 4000 years (consistent with the diameter of 
DEM L 71 and that the ejecta have all passed through the 
reverse shock as found by HGRS03), T cross /T ^>3. This 
suggests that during the current epoch at least, it is pos- 
sible for parts of the blast wave on opposite sides of the 
remnant to evolve independently in the manner we have 
assumed. 

10. LIMITS ON COSMIC RAY ACCELERATION 
PARAMETERS 

We can use our derived limits on the precursor thick- 
ness, preshock density, shock speed and explosion en- 



ergy to constrain cosmic ray acceleration parameters for 
DEM L 71 . These parameters are the cosmic ray dif- 
fusion coefficient, k, the maximum cosmic ray energy 
EreV and the roll-off frequency v ro u of the synchrotron 
emission spectrum. In the limit where k scales linearly 
with particle energy (the Bohm limit), the maximum cos- 
mic ray energy is E TeV « 1.1 (B /3fj,G) (E 51 V S /no) 1/3 
(Draine & McKee 1993), where B is the preshock mag- 
netic field strength and Vs is in kms -1 . Using our shock 
parameter estimates from the previous section, we find 
EreV < 11.3 (Bo/3 /uG). The roll-off frequency, defined 
as the frequency where the synchrotron power law emis- 
sion from cosmic rays begins to exponentially decline, is 
given by v roU = 1.4x 10 13 E^ eV (B /3 fiG) Hz (Reynolds 
1998). Inserting our estimate for Ex e v, we obtain v ro u < 
1.8xl0 15 (Bo/3 /zG) Hz, corresponding to a wavelength 
Koii > 1700 (3 fiG/B ) A. Furthermore, the cosmic ray 
diffusion coefficient is related to the precursor thickness 
via dp = k/Vs (Draine & McKee 1993). This gives an up- 
per limit k < 5.4xl0 25 //#+ cm 2 s _1 for cosmic rays in 
DEM L 71 . 

In an independent study of DEM L 71 , Hendrick & 
Reynolds (2001) extrapolated the observed radio syn- 
chrotron flux of DEM L 71 to the X-ray regime in order 
to estimate the maximum cosmic ray energy allowed by 
the data. They found Ex e v = 80^/10 fiG/B, where B is 
the postshock magnetic field strength. If we assume that 
B = 4 Bo = 12/iG, then we find that our estimated maxi- 
mum energy and roll-off frequency lie well within the limit 
allowed by the estimate of Hendrick & Reynolds (2001). 

In a study of 25 LMC supernova remnants by Mathew- 
son ct al. (1983), DEM L 71 was found to exhibit the 
lowest radio flux density in the sample (< 0.02 Jy at 408 
MHz) . This suggests that although DEM L 71 accelerates 
cosmic rays to TeV energies, the density of these particles 
is significantly lower than that found in the other LMC 
SNRs. Therefore, it seems unlikely that the blast wave 
structure of DEM L 71 is significantly altered by nonlinear 
feedback from cosmic ray acceleration. Further evidence 
against efficient cosmic ray acceleration in DEM L 71 may 
be found in our optical and Chandra data. The detailed 
correlation between the optical and X-ray morphologies 
(HGRS03, RGH03) suggests that the blast wave emission 
is predominantly thermal in origin. In addition, the shock 
models of Boulares & Cox (1988) predict precursor tem- 
peratures ~10 6 K for non-radiative shocks dominated by 
cosmic ray pressure. This is clearly ruled out by the Ha 
narrow component width measurements of SRL94 and the 
presence of Balmer line emission from behind the shock. 

11. SUMMARY AND DISCUSSION 

We have conducted a thorough spectral imaging study 
of the supernova remnant DEM L 71 in Ha. Spectra ex- 
tracted from the blast wave show the broad and narrow 
component emission lines characteristic of non-radiative 
shocks in partially neutral gas. The Ha broad compo- 
nent width ranges from ~ 450 km s _1 along the eastern 
edge to values ~ 1000 kms -1 along the faint western edge. 
No broad component is detected in spectra of the faint 
northern rim, indicating a very low density in these regions 
and probably a broad component width in excess of 1000 
kms -1 . The range of shock velocities is 500 kms -1 along 
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the eastern rim, 700—800 km s" 1 along the southeastern 
and southern rims and 800—1000 kms -1 along the west- 
ern and northwestern rims, other parts of the blast wave. 
The Ha broad-to-narrow ratios Ib /In range from around 
0.4 to 0.8, falling significantly below the values predicted 
by our shock models. The most likely explanation for the 
discrepancy is extra narrow component emission from a 
cosmic ray precursor, although a fast neutral precursor 
may also contribute some emission. 

If collisional excitation in a 40,000 K cosmic ray pre- 
cursor produces narrow line emission, the steep Balmer 
decrement of the precursor should affect the H/3 broad-to- 
narrow ratio far less than the Ha broad-to-narrow ratio. 
Therefore, a measurement of the H/3 broad-to-narrow ratio 
in the future would be highly desirable. Among the other 
studied Balmer-dominated remnants, Tycho also exhibits 
an anomalously low Ha broad-to-narrow ratio (SKBW91, 
Ghavamian 1999, Ghavamian ct al. 2001). In that case, 
where both Ha and H/3 broad-to-narrow ratios were mea- 
sured, the former ratios were found to lie 40% below the 
lowest predicted values, while the latter ratios were found 
to agree with the lowest predicted values. This is consis- 
tent with detection of superimposed precursor and shock 
emission. Note that the high inferred preshock neutral 
fraction for Tycho's SNR (GRHB00) would make the pre- 
cursor emission more detectable there than remnants like 
SN 1006, where the inferred preshock neutral fraction is 
only ^10% (Ghavamian ct al. 2002). A high preshock 
neutral fraction and large size for DEM L 71 (over twice 
that of Tycho) would certainly elevate the role of a pre- 
cursor in understanding the blast wave emission. 

Our ability to map the shock velocity along the entire 
rim of DEM L 71 has enabled us to estimate the age and 
explosion energy of this remnant entirely from the blast 
wave kinematics. Assuming the blast wave propagates into 
a radially uniform ISM, the derived age for DEM L 71 is 
3500—5000 years, less than half the value calculated from 
previous observations (SKBW91). The explosion kinetic 
energy is (2 — 6) x 10 50 ergs, less than half the canonical 
value assumed for Type la SNe. One possible explana- 
tion for the difference may be the loss of energy during an 
earlier phase of evolution, when DEM L 71 experienced 
efficient cosmic ray acceleration. However, the lack of ev- 
idence for nonthermal emission in Chandra X-ray spectra 
of the blast wave (RGH03), the strong similarity between 
the optical and X-ray morphologies of the blast wave and 
the faintness of the radio emission (Mills 1984) suggest 
that blast wave energy loss to cosmic ray acceleration is 
not significant in the current epoch. 

The presence of significantly neutral gas around the en- 
tire blast wave of DEM L 71 places a strong limit on the 
amount of photoionizing radiation emitted by the accret- 
ing white dwarf progenitor. SuperSoft X-ray Sources (SSS) 



are believed to be one class of these progenitors (van den 
Heuvel et al. 1992, Rappaport, Di Stefano & Smith 1994). 
In their detailed calculations of ionization conditions sur- 
rounding SSS, Rappaport et al. (1994) found that H is 
entirely ionized inside a radius of 28 pc for a source with 
standard parameters T e ff = 4x 10 5 K, L = 10 38 ergs s _1 in 
a medium with n = 1 cm~ 3 . This is clearly larger than the 
current blast wave radius of DEM L 71 (R < 9.5 pc) and 
limits the luminosity of an SSS progenitor to values < 10 37 
ergs s _1 . Note that due to the low ambient density, the 
recombination time per H atom is ~ 10 5 yr, so that lit- 
tle neutral gas should re-form within the ionized bubble 
created by an SSS progenitor. Moreover, it is worth not- 
ing that SSS progenitors with L ;> 10 37 ergs sm _1 also 
appeared to be ruled out in younger Balmer-dominated 
remnants of Type la SNe, such as SN 1006 and Tycho's 
SNR. This limit is especially stringent in the latter rem- 
nant, where the blast wave encounters significantly neutral 
gas only 2—3 pc from the center of the explosion. 

There are at least 4 other SNRs with nearly com- 
plete Balmer-dominated rims, namely RCW 86 (Long & 
Blair 1990, Smith 1997) in our galaxy and 0509-67.5, 
0519-69.0 and 0548-70.4 (Tuohy et al. 1982) in the 
LMC. Of these remnants, 0509—67.5 shows no evidence 
of a broad component (Tuohy et al. 1982, SKBW91), 
probably due to its very high shock speed. On the 
other hand, although spectra of RCW 86 have revealed 
strong broad component emission (Long & Blair 1990, 
Ghavamian et al. 2001), it is nearly 40' in size (Smith 
1997). This is larger than the FOV of most Fabry- Perot 
spectrometers, making a detailed FP study of this SNR 
more time consuming and challenging. Only 0519—69.0 
and 0548—70.4 appear to be small enough (< 1') with 
small enough broad component widths (750-1300 kms -1 , 
SKBW91) to be amenable to the type of study we have 
undertaken of DEM L 71 . By applying the blast wave 
kinematic analysis described in this paper to 0519—69.0 
and 0548—70.4, we would be able to enlarge our sample of 
SNRs with well determined ages and explosion energies. 
Such a study would be invaluable in testing evolutionary 
models of young SNRs and in the interpretation of existing 
data from the infrared to the X-ray regimes. 
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Fig. 1. — The locations, sizes and orientations of spectral extraction apertures for the DEM L 71 blast wave. The image here is centered on 
the rest velocity of the narrow component Ha line. The width of each extraction aperture is 3''2. The two black boxes near the center mark 
the lcr error bars on the derived geometric center from fits to Vs(R) (larger box) and nrj(-R) (smaller box). The overlapping region between 
the two is 0"7 X 2"6 in size. For comparison, the expansion center of the iron-rich ejecta estimated from the Chandra data (HGRS03) is 
marked by the cross. North is at the top, east is to the left. 
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Fig. 2. — The sky-subtracted Ha line profiles of blast wave apertures f through 4. The data are marked by filled circles with error bars. 
The total line profile fits are shown along with the individual broad and narrow profiles. The Ha lines from apertures f and 2 are well fit by 
a single narrow component profile, suggesting that the broad component is either nonexistent or too broad and faint to be detected. 
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Fig. 3. — The sky-subtracted Ha line profiles of blast wave apertures 5 through 8. The total line profile fits are shown along with the 
individual broad and narrow profiles. 
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Fig. 4. — The sky-subtracted Ha line profiles of blast wave apertures 9 through 12. The total line profile fits are shown along with the 
individual broad and narrow profiles. 
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Fig. 5. — The sky-subtracted Ha line profiles of blast wave apertures 13 through 16. The total line profile fits are shown along with the 
individual broad and narrow profiles. 
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Fig. 6. — The broad-to-narrow ratios predicted by shock models for combined blast wave aperture regions. The shaded regions indicate the 
observed range of Ha broad-to-narrow ratios. The curves are computed for a range of (T e /T p )o at the shock front, for two selected preshock 
H neutral fractions: 50% neutral H (dot-dashed curve) and 90% neutral H (solid curve). 
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Fig. 7. — Correlation between the blast wave speed on radius in DEM L 71 , shown for the combined apertures. The best linear fit to the 
Scdov relation Vg = 2 11 ; s marked by the solid line. The la limits on the slope are indicated by the dashed lines. 
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Blast Wave Radius (pc) 

Fig. 8. — Correlation between the total preshock density and radius for DEM L 71 . Preshock densities for the combined apertures are 
determined from Chandra spectral models of RGH03. The data are fit to the Sedov relation no(R) oc R~ 5 (sec Equation 5 in the text). The 
solid line shows the best fit curve, while the dotted lines denote the la limits on the fit. 
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Table 1 

Blast Wave Extraction Apertures 



Aperture 


a(2000) 


(5(2000) 


Size (") 


Position Angle (o) a 


1 


05:05:42.0 


-67:51:57.5 


3.2x12.3 


69 


2 


05:05:40.2 


-67:52:02.6 


3.2x11.0 


69 


3 


05:05:38.6 


-67:52:08.5 


3.2x9.7 


55 


4 


05:05:37.7 


-67:52:16.2 


3.2x8.4 


30 


5 


05:05:35.8 


-67:52:42.8 


3.2x9.7 





6 


05-05-35 8 


—67-52-53 8 


3 2x12 3 


o 


7 


05:05:37.7 


-67:53:10.1 


3.2x12.3 


115 


8 


05:05:39.5 


-67:53:14.4 


3.2x9.7 


115 


9 


05:05:41.2 


-67:53:18.5 


3.2x11.0 


115 


10 


05:05:43.2 


-67:53:19.7 


3.2x11.0 


90 


11 


05:05:45.5 


-67:53:13.5 


3.2x11.0 


35 


12 


05:05:46.4 


-67:53:05.4 


3.2x8.4 


35 


13 


05:05:46.9 


-67:52:30.1 


3.2x7.1 


340 


14 


05:05:46.3 


-67:52:22.9 


3.2x8.4 


330 


15 


05:05:37.9 


-67:52:55.9 


3.2x9.7 


320 


16 


05:05:41.4 


-67:53:09.1 


3.2x11.0 


110 



a Measured E of N 



Table 2 

Fitted Parameters for Blast Wave Ha Profiles 



Aperture 


X 2 / 9 d.oi. a 


P(>F) b 


Vfwhm (kms 1 ) c 


Ib/In 


AV B - N d 


1 


1.17 










2 


0.56 










3 


0.92 


0.05 




0.50±0.08 


+851°? 


4 


2.59 


0.05 


700+^ 


0.51±0.06 


+27±» 


5 


1.10 


0.04 


1055118 


0.88l°;lJ 


-16±g 


6 


0.48 


0.10 


9oo;^ 5 5 


0.20±0.08 




7 


0.54 


0.07 


59511°° 


0.29±0.09 


-591^° 


8 


1.26 


0.05 


83011^ 


^+°- 09 

u - oo -0.08 


,o+50 
~l~ z -54 


9 


1.09 


0.05 


83511^ 


r;i -4-0.09 
u - ol -0.08 


-103±H 


10 


1.16 


0.11 


3501^° 


0.26±0.09 


+104+™ 


11 


1.07 


0.07 


oyo -100 


0.36±0.09 




12 


1.17 


0.04 


7851*5 


9S+ 011 


-74lf 7 


13 


0.60 


0.04 


5401*° 


46+°- 07 


+21^ 


14 


2.53 


0.08 


355±90 


4^+ - 16 


i o+23 
io -28 


15 


0.81 


0.06 


4901^ 


0.28±0.07 


+160l°| 


16 


0.31 


0.04 


86511*° 


0.67±0.11 





a The quoted % 2 values are for one-component profile fits for apertures 1 and 2 and two-component profile fits for 
apertures 3-16. 

b P(>F) denotes the F-test goodness-of-fit for inclusion of a broad component (see text for details). 
c The Gaussian FWHM of the broad component. 

d The shift of the broad component centroid from the narrow component centroid. 
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Table 3 

Ha Profile Fits For Combined Blast Wave Apertures 



Aperture 


X 2 /9 d.o.f. 


P(>F) a 


Vfwhm (kms 1 ) b 


Ib/In 


3+4 


1.87 


0.05 




0.51+0.06 


5+6 


0.68 


0.04 




0.54+0.09 


7+8+9 


0.73 


0.04 


80511S 


n 4Q+ - 07 


11+12 


0.40 


0.04 


735lg° 


0.66+0.08 


13+14 


1.90 


0.04 


450+60 


n 4 4 +0.06 
u - 44 -0.05 



a P(>F) denotes the F-test goodness-of-fit for inclusion of a broad component (see text for details). 
b Thc Gaussian FWHM of the broad component. 



Table 4 

Derived Shock Velocities for Combined Blast Wave Apertures 







(%)o = l 




Chandra Estimates* 




Aperture 


V s (km s" 1 ) 


V s (km s- 1 ) 


V s (km s- 1 ) 


n (cm -3 ) 


(5f)o 


3+4 (Xl) b 


01 K +115 
010 -100 


1055 ±iS8 


855 US 


n 4t;+0-08 

U - 40 -0.12 


24 +o.i2 


5+6 (X2) 


960±?« 


12401*?° 


980^225 


n oo+0. 05 
U - JJ -0.10 


0.1+0.29 


7+8+9 (X3) 


77c +140 
' 10 -115 


1005 


77 c +145 
' '°-110 


o-48±g-.2g 


0.01+003 


11+12 (X4) 


7io±Sg° 


QIC +130 
yl ° -105 




0A7t° o f 5 


o.4fi±g:ig 


13+14 (X5) 


430 t«j 


555 7q 


555 


1.30±° ; 2 2 ° 


1.0-0.16 



a (T e /T p )o is derived from X-ray spectral models of RGHW02. The appropriate shock speed for that equilibration is quoted in column 4. 
'Names of the corresponding X-ray spectral extraction regions used in the Chandra analysis of RGHW02. 



